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ABSTRACT 

We consider the nature of orbits near the solar neighborhood which are perturbed by local spiral arms 
and the Milky Way bar. We present a simplified Hamiltonian model which includes resonant terms from 
both types of perturbations and is similar to the forced pendulum. Via numerical integration of this 
model we construct Poincare maps to illustrate the nature and stability of the phase space. We find that 
resonance overlap is most likely to cause widespread chaos when the pattern of the spiral structure puts 
the solar neighborhood near the 2:1 inner Lindblad resonance (ILR) in the case of a 2-armed pattern, 
or near the 4:1 ILR in the case of a 4-armed pattern. When this happens both the quasiperiodic orbits 
which support the spiral structure and those that oscillate with the bar are disrupted near the bar's 
2:1 outer Lindblad resonance (OLR). Consequently the pattern speed of spiral structure which passes 
through the OLR must be faster than ~ 0.45 times the solar neighborhood angular rotation rate if it is 
2-armed or faster than 0.75 times this value if it is 4-armed. Alternatively the OLR may form a boundary 
between spiral modes at different pattern speeds. In all cases we find that spiral structure is disrupted 
by the OLR over a narrow range of radius and the extent of the orbits aligned perpendicular to the bar 
at the OLR is limited by the spiral perturbations. We find that the boundary at an asymmetric drift 
velocity of v ~ —30 km/s of the w-anomaly in the velocity distribution in the solar neighborhood is due 
to the abrupt bifurcation of the orbit families associated with the OLR. The upper boundary at v ~ —60 
km/s is truncated by the spiral structure. The radial velocity width of the anomaly is probably bounded 
by chaotic regions which surround the family of quasiperiodic orbits oriented perpendicular to the bar. 

Subject headings: 



I. INTRODUCTION 

It has been established, beyond doubt, that the Milky 
is barred as are many nearby galaxies. Evidence for the 
bar comes from asymmetry in the bulge surface photom- 
etry, star counts and interpretation of the gas kinematics 
(for a recent review see Gerhard 2002). The most common 
values suggested from these observations are a corotation 
radius of 3.5-5 kpc and a major axis in-plane angle with 
respect to the Galactic Center of 4>bar ^ 15 — 45°. 

The solar neighborhood local velocity distribution is ex- 
pected to contain structure caused by large scale nearby 
gravitational perturbations. While structure at low asym- 
metric drift is likely to be affected by spiral structure and 
disrupted stellar clusters (e.g., Skuljan et al. 1999), the 
Hercules stream at a mean heliocentric asymmetric drift 
velocity of v ~ 45km/s is suspected to be related to the 
2:1 outer Lindblad resonance (OLR) with the Galactic bar 
(Dehnen 1999; Fux 2001; Raboud et al. 1998). This stream 
is also known as the it-anomaly. Due to the older and 
rather high metallicity stars in it, the stream probably 
has a stable kinematic origin (Raboud et al. 1998; Dehnen 
1998). The local velocity distribution as constructed from 
Hipparcos catalog most clearly shows this stream as a 
strong and separate feature Dehnen (1998); in fact there 
is a saddle point in the distribution between the stream 
and the bulk of the stars at v ~ 30 km/s (see Figures in 
Dehnen 1998; Fux 2001). 

Linear theory predicts that the orientation of orbits will 
shift across the OLR, from oriented along the bar major 



axis outside the OLR to perpendicular to the bar within 
it (Binney & Tremaine 1987). Near the peak of the res- 
onance both types of orbits can exist (Contopoulos 1975; 
Weinberg 1994). Dehnen (2000) showed using a backwards 
integrating bar growth model that stars were likely to be 
captured into resonant islands associated with these 2 ma- 
jor orbit families. However the final velocity distribution 
is strongly dependent upon the timescale over which the 
bar grows as well as the assumed initial velocity distri- 
bution. Fux (2001) considered the stability of orbits and 
suggested that chaotic regions associated with resonances 
were likely to cause over-densities in some regions of phase 
space and under-densities in others. He illustrated with 
N-body simulations that the resulting stellar distribution 
could correspond to that observed. 

In some sense it is surprising that the OLR with the 
bar provides such a good explanation for the u-anomaly 
because locally both gaseous and stellar tracers primar- 
ily show spiral structure near the solar circle (e.g., Vallee 
1995; Drimmel & Spergel 2001). Previous works have not 
considered the more complicated problem caused by the 
forces due to local spiral structure in addition to those 
caused by the bar. 

In this paper we consider the dynamics of stars that are 
affected by perturbations from both spiral structure and 
the Milky Way bar. We construct a simple Hamiltonian 
model for the dynamics which contains resonant terms at 
both the bar pattern speed and a slower pattern speed 
from local spiral structure. We find that the Hamilto- 



1 



2 



nian can exhibit resonance overlap and so widescale chaos. 
Consequently we address the issue of orbit stability by 
computing area preserving or Poincare maps for different 
values of the spiral pattern speed. 

2. HAMILTONIAN FORMALISM FOR THE KINEMATICS 

As shown by Contopoulos (1975, 1988), the dynamics 
of stars confined to the galactic plane moving in a smooth 
Galactic potential lacking non-axisymmetric perturbations 
is described by a Hamiltonian which can be written in a 
third order post epicyclic approximation as 

H (I 1 ,e 1 ;I 2 ,0 2 ) = h + Luxh + uj 2 h + al\ + 2bIJ 2 + cl\... 

(1) 

The high order of the expansion is required to exhibit the 
bifurcation of the major resonances. The action variables 
I\ = j- / rdr and I 2 = J — J c are integrals of motion 
when the Hamiltonian is unperturbed. Jo is the particle's 
angular momentum and J c the angular momentum of a 
particle in a circular orbit at a radius r c which is the ra- 
dius of a circular orbit with energy h. The frequencies 
and constants in the above Hamiltonian are evaluated at 
this radius. We can either to choose to work in a frame 
rotating with a perturbation pattern or an inertial one. 
Contopoulos (1975) worked in a frame rotating at the bar 

or spiral pattern speed, Q s , so h = + Vo(r c ) — Q s J c and 

oj 2 = fi c — fig where f2 c is the angular rotation rate of a 
circular orbit. Since we will consider more than one pertur- 
bation and each perturbation will have a different pattern 
speed, in this paper we choose to work in an inertial frame. 

Consequently h = t^St + Vb(r c ) and uj 2 = ^c- In either 

frame, J c = r^,Q, c is the angular momentum and Voir) the 
axisymmetric component of the potential. The frequency 
u>\ = k c where k c is the epicyclic frequency which is eval- 
uated at r c . Unfortunately the value of r c depends on Q s 
so transferring coordinate systems between an inertial and 
a rotating one is not trivial. 

The constants a, b, c are given in the appendix of Con- 
topoulos (1975). When the rotation curve is flat k c = 
\An c , a = b = and c ; 
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The radius r, the angle 9 and their time derivatives are 

related to the action angle variables to first order in If 
(from the appendix of Contopoulos 1975) 
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where Q\ is the epicyclic angle and 9 2 is the azimuth of the 
epicyclic center. 6\ and 9 2 are the angle variables conju- 
gate to I\ and I 2 . 

To consider the affect of non-axisymmetric perturba- 
tions in the gravitational potential caused by a bar or spi- 
ral structure we must estimate the strength of these per- 
turbations in terms of the action angle variables described 
above. We concentrate on the Lindblad resonances for a 



number of reasons. They are first order in the amplitude 
of the perturbation and in I\ 5 and so likely to be strong. 
The 2:1 OLR from the Galactic bar is near the location of 
the solar neighborhood. Stars on orbits influenced by the 
ILRs (inner Lindblad resonances) are required to self con- 
sistently support spiral structure (e.g., Contopoulos 1988). 

2.1. The perturbation from the Milky Way Bar 

For a bar like perturbation we assume that the non- 
axisymmetric component of the gravitational potential de- 
pends on radius and can be expanded in Fourier compo- 

I1C lit s 

V 1 {r,9) = Y,B m (r)cos[m{0-n b t)} (6) 

m 

where fib is the angular rotation rate of the bar. In action 
angle variables we can write the potential perturbations 

1 12 

to first order in I 1 and the strength of the perturbation 
(following expressions given by Contopoulos 1988) 
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where the first and second angular terms correspond to the 
m : 1 outer and inner Lindblad resonances respectively. 
Here 

(i m = T. B ' m + ~r~ Bm ( 8 ) 

_ ^ ' c™c 

where B m and B' m are evaluated at r c . 

The solar neighborhood is well past the end of the Galac- 
tic bar so we can approximate the potential perturbation 
as a quadrupolar term in the gravitational potential. 
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where r& is the radius at which the bar ends, r& w 0.45ro, as 
found from IR photometry (e.g., Dwek et al. 1995) where 
r is the radius of the solar circle from the Galactic Cen- 
ter. The pattern speed of the bar is constrained by to 
be about 1.85 times the solar neighborhood value of the 
angular rotation rate, flo (Dehnen 1999). Fux (2001) con- 
sidered bars about twice as strong as those considered by 
Dehnen (1999). 

Dehnen (2000) estimates a = ^ (n^j « 0.01, so 
at the solar circle ^ <~ 0.036. From equation (8) and 
assuming that the rotation curve is flat we find that 

j3 2 ~ 0.086^ (jr^j so at the solar neighborhood we es- 
timate 



02 ^ 0.003-^ (10) 
where v c is the velocity of a star in a circular orbit. 

2.2. The perturbation from local spiral structure 

For spiral structure we assume that the radial variations 
depend on angle and the amplitude is nearly constant with 
radius. 
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Here to refers to the number of arms, fl s to the spiral pat- 
tern speed, k m to the wavenumber and a m to an angular 
offset. 

1 /2 

Contopoulos (1975) showed to first order in I 1 ' the po- 
tential perturbation 

V hm (I u e 1 ;l2,e 2 ) = e m x 

[cos(6»i - to(0 2 - O s t + 7 m+ )) + cos(0i + to(0 2 - fi s t + 7 m 
where 



(13) 



for k J n <C 1. The angular offsets j m ± depend on a m and 
the wavenumber (see Contopoulos 1975 for expressions). 

Tracers of spiral structure in the solar neighborhood sug- 
gest that the spiral structure is quite tightly wound (e.g., 
Vallee 1995). In the WKB approximation 

A ~ - 2G ^ A (14) 

\r"m I 

where S m is the amplitude of the spiral surface density 
variations divided by the mean surface density (Binney 
& Tremaine 1987). For a tightly wound structure we ex- 
pect that e m <~ Am ^ m (taking the largest term in equation 
13). The mean surface density of disk mass in the so- 
lar solar neighborhood is S ~ 50M Q pc~ 2 (Holmbcrg & 
Flynn 2000), and we use v c ~ 200km/s. A two armed spi- 
ral is seen in the near infrared COBE/DIRBE data with 
pitch angle in the range p ~ 15.5° - 19° (Drimmel & 
Spcrgcl 2001), giving k m ro = mcotp in the range 5.8 - 
7.2. The amplitude of the spiral structure could be as 
large as S m - 0.5 - 1 (Drimmel & Spergel 2001). We 
estimate for the solar neighborhood 



obtaining a resonant angle and new momenta 

</> = e 1 + 2(e 2 - n b t) 

h = 2J 1 + J 2 . (18) 
Ji = Ji and 9 2 remain unchanged. Our Hamiltonian be- 
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The expression for e m is independent of k m and primarily 
depends on the amplitude of the spiral density variation. 

3. HAMILTONIAN WITH PERTURBATIONS FROM BOTH 
SPIRAL STRUCTURE AND THE MILKY WAY BAR 

We now take the perturbations due to the bar and spiral 
arms and add them to the unperturbed Hamiltonian. We 
have estimated the size of the coefficients (3 and e for these 
perturbations in the previous section. 

We expect that the 2:1 OLR with the Galactic bar and 
the 2:1 or 4:1 ILR with local spiral structure will be the 
strongest resonances near the solar neighborhood. So ig- 
noring all other resonant terms we can simplify the prob- 
lem to 

H = all + kIi + VI 2 + ft^)* cos[0i + 2(0 2 - n h t)] 

+e m (2ii)s cos[0i - m(9 2 - fl s t) - 7] + .... (16) 

where to = 2 for the 2:1 ILR and to = 4 for the 4:1 ILR. 
We have assumed that I 2 is small and so have dropped the 
terms bl\l 2 and c/| from the Hamiltonian. 

We do a canonical transformation with generating func- 
tion 

F 2 (j u j 2 ,e u e 2 ) = (0i + 2(0 2 - n 6 t)) + j 2 9 2 (17) 



where 



5 = K + 2(n-n b ) 

a = (mQ s + 2Vl b )t - (to + 2)0 2 . 

We can approximate 02 f=s fit so that a w vt where 

v = mfl s + 2tt b - (2 + m)fl. 

We rewrite our Hamiltonian in a simpler form 

H(J U <P) w aJ\ + SJ! - cos(^) 

-e 2 (^-)5cos((/> + ^-7). 

We rescale by dividing by a, put lengths in units of r c and 
time in units of i so that J\ = jiKr 2 . Since a < this has 
the effect of reversing the sign of t. 

„ _ _ 1 I 

Kh,4>) ~ Ji+Sji-(3 2 ji cos(0)-e 2 Ji cos(</>-^-7) (24) 
where the unitless coefficients 

k + 2(n - n b ) 
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(25) 
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If we set e = then h = jf+Sji —f3j^ cos <f>. This Hamil- 
tonian is in the same form as e — e resonances discussed 
in the context of solar system orbital resonances (Murray 
& Dcrmott 1999; Borderies & Goldreich 1984). As illus- 
trated by Murray & Dermott (1999); Borderies & Goldre- 
ich (1984) the resonance bifurcates at a critical value of a 
parameter which depends on 5 and (3 (see Figure 1). The 
resonance has three fixed points when 

i < "I- (29) 
3pa 

One of the fixed points is unstable, the other two cor- 
respond to two resonant stable islands of libration at 
4> = 0, 7r. Only one fixed point exists when the above 
inequality does not hold, and its location is determined by 
the sign of (3. When (3 > this fixed point is at <j> = 0, 
other wise _at <j> = ir. 

When j3 = the Hamiltonian can be put in the form 
h = jf + (5 — D)ji — e?2 coscp' with a suitable canonical 
transformation. In this case the bifurcation happens when 



2(5-9) 
3ei 



< -1. 



(30) 



The distinction between regions with two stable fixed 
points and with one fixed point is important when reso- 
nances overlap. For this type of Hamiltonian a separatrix 
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exists only when there are two stable fixed points. Be- 
cause the period of the orbits becomes infinite at the sep- 
aratrix, additional perturbations are most likely to cause 
instability and a chaotic region near the original location 
of the separatrix. When we have two resonant terms in 
the Hamiltonian and the resonances overlap, large regions 
of phase space will become strongly chaotic only when at 
least one of the resonances contains a separatrix. 
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Fig. 1. — We plot level contours of the Hamiltonian 
4>) = j 2 + jS — (Sj 1 / 2 cos <f>. Dashed contours are 
negative. The axes are x,y are defined by Eqn(31). 
The critical value for the resonance to bifurcate hap- 
pens in the middle panel where = — 1. The top 

3/33 

of the panel contains a fixed point which would cor- 
respond to a closed or periodic orbit aligned with the 
bar so that <j> = 0. This is the situation outside the 
OLR. Inside the OLR the resonance bifurcates, and 
both periodic orbit families are present. This situa- 
tion is shown in the bottom panel. The fixed points 
at (f> = n correspond to periodic orbits oriented per- 
pendicular to the bar. Only when < — 1 is there 

3/33 

a separatrix. Additional perturbations are most likely 
to cause instability when there is a separatrix. 

We now consider the orientations of the orbits and signs 
of the individual terms. Fixed points in our simple one di- 
mensional system correspond to closed or periodic orbits 
in the full 2 dimensional system. Assuming that Qb is 1.85 
times the angular rotation rate at the solar circle (Dehncn 
1999) and a flat rotation curve we estimate 5 ~ 0.2 at the 
solar circle. This drops to zero as we approach the 2:1 OLR 



at a Galactic centric radius of r w 0.9ro- Because a ~ — 1, 
§ ~ fo^c/Vc an d e ~ er c/vc- For radii outside the OLR, 
8 > and we expect only one fixed point or periodic orbit. 
Because > the quasipcriodic orbits at the solar circle, 
outside the OLR, will be aligned with the bar. They are 
referred to as the Xi(l) orbits (e.g., see Fux 2001; Dehncn 
2000). Inside the OLR the resonance bifurcates and two 
families of periodic orbits (fixed points) exist, the Xi(l) 
and xi(2) families. The xi(2) family is aligned perpen- 
dicular to the bar and found at <p = 7r. We expect large 
scale chaos can occur when there is a separatrix, 5 < 0, 
and inside the radius of the OLR. 

fis is ~ 0.3 times the angular rotation rate at the solar 
circle when r is just outside the 2:1 ILR and 0.6 if r is 
just within the 4:1 ILR. 7 is related to the phase of the po- 
tential for the spiral pattern, and is difficult to constrain 
from observations. For spiral structure at our location in 
the solar neighborhood we should be outside the m:l IRL 
of the pattern so that 5 — v < 0. If we aren't outside 
the resonance, the stellar orbits will not support the spi- 
ral structure, and will negate it instead (e.g., Contopoulos 
1988). We expect e < for periodic orbits with = 7 
to be in phase with the spiral arm. By in phase we mean 
that the orbit will be aligned with radial maxima on the 
same axis as the density maxima. 

3.1. Generating Poincare maps 

Previous works addressing the stability of the orbits in 
the solar neighborhood have computed Liapunov expo- 
nents, studied N-body simulations (Fux 2001) or carried 
out backwards integration schemes (Dehnen 2000). Here 
we adopt the approach of Fux (2001) and strive to identify 
regions which regions of phase space that can keep stars 
for long periods of time. However as we explain below we 
do this by mapping the nature of phase space instead of 
computing Liapunov exponents. 

The Hamiltonian (Eqn.24) contains two main resonant 
terms separated by a frequency, v. This is similar to the 
forced pendulum. When the two resonances overlap, a 
chaotic region can be generated at the separatrix of one of 
the resonances. This zone has a Liapunov time ~ (Hol- 
man & Murray 1996). Outside the main chaotic zone we 
expect quasiperiodic islands. This picture is qualitatively 
different than that explored by Fux (2001) who illustrated 
narrow bands of unstable regions, each with different Lia- 
punov timescales. We expect instead large bands of chaos 
described by one Liapunov time which arc surrounding 
stable islands. Consequently we do not compute Liapunov 
times for a range of initial conditions but instead map the 
structure of the phase space itself. 

To investigate the stability of the system containing per- 
turbations from both the bar and spiral structure we inte- 
grate Eqn(24) numerically. This system is time dependent 
so the Hamiltonian itself is not conserved. However if we 
plot points every timestep P„ = ^zl, we derive an area pre- 
serving or Poincare map. This procedure generates maps 
which are like surfaces of section, and these we can use 
to study the stability of the orbits. In such a map orbits 
are cither area filling or curved loops. We denote area 
filling orbits as chaotic (in the sense that orbits diverge 
exponentially) and the curved orbits as quasiperiodic. 

This procedure is excellent for addressing the question 
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of orbit stability and identifying the regions accessible to 
individual orbits. But because we plot every P„ we can- 
not tell if an orbit is librating around <j> or around an angle 
<fi — vt. This makes it difficult to determine whether fixed 
points are associated with (or supporting) the bar or the 
spiral pattern. However we can qualitatively deal with this 
problem by changing the orientation of the spiral pertur- 
bation. We do this by adjusting 7. We set 7 = | so that 
orbits in phase with the spiral pattern are located at an 
angle of <j> ~ ^. Orbits perpendicular to the spiral pat- 
tern are located at an angle of — | . In comparison, orbits 
associated with the bar's OLR librate about (p = or n. 

As is commonly done, (e.g., Murray & Dermott 1999; 
Bordcries & Goldrcich 1984), we plot all figures in this 
paper in a coordinate system with 
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(31) 



With respect to our action variable, j\ — s 2 /2 where s 
is the radial distance on these plots. Since we defined 
•h = ji Kr l, the cpicyclic amplitude is approximately the 
same as s, the radial distance in these plots. In other 
words s gives the maximum difference between the radial 
position of a particle and r c . 

For each series of integrations (Figures 2-4) we assume 
a value for f2 s and fif, (in units of the_solar neighborhood 
angular rotation rate flo, and e and j3 which are directly 
estimated from the strengths of the spiral structure and 
bar (see previous sections). We define the radial offset, 
dr, from the radius of the sun as 



dr = 



ro 



(32) 



so that = ^ = 1 + dr. S and v for each integration are 
calculated from O s and fit assuming a flat rotation curve 
and for a range of dr. In the figures the different panels 
correspond to different values of dr. 

Considering the bar only, as we vary dr and so 5, phase 
space changes in appearance from that shown at the top 
of Figure 1 to that in the bottom of this figure as we pass 
through the OLR. Considering the spiral arms only, phase 
space should look like the bottom of Figure 1 but rotated 
by 90 degrees since we choose 7 = it/2. When S < 0, 
the bar OLR gains a separatrix which is likely to become 
unstable when perturbed by the spiral arms. 

In Figure 2 we show the result of integrating Eqn.(24) 
for a 2-armed pattern near the 4:1 ILR of the spiral pattern 
for moderate spiral and bar strength. For each of 30 parti- 
cles chosen with different initial conditions, 200 timesteps 
are plotted with dt = 2n/v. For the panel on the top left, 
corresponding to the Galactic radius of the sun, the orbits 
librate about <j> = w/2 and support the spiral structure. At 
dr <~ —0.90, orbits become are aligned with the bar (xi(l) 
type orbits) and librate about </> = 0. At dr = —0.095 the 
xi(2) orbits appear. These librate about <fr = — 7r and are 
oriented perpendicular to the bar. As dr decreases fur- 
ther the orbits librate closer and closer to <j> = tt/2 and 
so support the spiral structure and are no longer aligned 
perpendicular to the bar. For the intermediate region be- 
tween dr = —0.10 and —0.175 the quasiperiodic orbits are 
likely to be oscillating with both the bar and the spiral 
structure. 



Just past the OLR (the center of the resonance is near 
6 = 0, but a separatrix appears only at negative values for 
(5), we see in Figure 2 that the quasiperiodic regions are 
bounded by thick bands of area filling or chaotic orbits. 
The spiral structure is disrupted (unsupported) between 
dr = -0.075 and -0.115. Past dr = -0.175 the Xl (2) 
orbits are limited by the influence of the spiral structure. 
In short the spiral structure limits the extent of the orbits 
perpendicular to the bar OLR and the bar OLR disrupts 
the spiral structure over a narrow range of radius. 

3.2. On the u,v plane 

To sec if our maps correspond to what is observed in the 
local velocity distribution we must consider what values of 
our action angle variables correspond to the u, v velocities 
that arc used to measure stars in the solar neighborhood. 

The solar neighborhood velocity distribution is typically 
plotted with respect to the asymmetric drift velocity (az- 
imuthal), v, and the radial velocity, u, where u > cor- 
responds to velocities toward the Galactic center. Orbits 
at lower v values have lower r c values and are expected to 
oscillate about galactic radii that are smaller than ro, the 
radius of the solar circle. We now estimate what values of 
dr correspond to the center and boundary of the Hercules 
stream or u-anomaly. 

If we assume a flat rotation curve, the Hamiltonian lack- 
ing non-axisymmetric perturbations 
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In r/r . 



n 2«2 2v 2 c 
For a circular orbit of radius r c both u, v = and 
H 



— = 1.0 + In r c /r . 



vf. 



(33) 



(34) 



We set r = ro at the solar neighborhood, equate the two 
expressions and solve for r c as a function of the asymmet- 
ric drift velocity v, finding In ^ = ^- + u . We expand 
In — and obtain 

ro 



dr = -(1 + _) + —. 

v c 2v c vt 



(35) 



Assuming u = 35km/s at the center of the stream and 
v c = 200km/s 



dr 
dr 



-0.14 
-0.18 



for v 
for v 



-30kr 
-45kr 



(36) 



where v <~ —45 km/s at the center of the Hercules stream 
or u-anomaly and v <~ —30 km/s at the edge or boundary 
of the stream. 

The stream itself is quite wide along the u direction 
(over 100 km/s wide) and comparatively quite narrow in 
the v direction (40 km/s wide). In fact the mean u value 
of the stream is much less than its u width. If we use a 
width in u of 35 ± 50km/s then we estimate a range for 
-0.20 < dr < -0.07 in the stream. 

In Figure 2 we see the x\(2) OLR orbits appear at 
dr ~ 0.1 and not at dr <~ —0.15 which is the location 
of the Hercules stream or w-anomaly. We could adjust the 
pattern speed of the bar to move the resonance over, in 
fact with il s = 1.95, we can match the location of the 
■u-anomaly. However before we do this we should con- 
sider the approximations made in our analysis. In the full 
Hamiltonian (Eqn.l) there is a term proportional to I1I2 
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which, had we kept it in our analysis, would have caused 
an additional term 7-7/72 added to S in Eqn.(26). In drop- 

CI K C T C 

ping the extra terms in the Hamiltonian we have assumed 
that I2 is small. Because b is the same sign sign as a, 5 
should be larger than we have calculated. This would have 
the affect of shifting the location of the resonance to more 
negative dr or larger distances from the solar neighbor- 
hood. 

The maps shown in Figures 2-4 do show the likely mor- 
phology of phase space near the OLR. However we should 
not use them to constrain the bar pattern speed at a level 
better than 25% unless we have taken into account this 
correction. We estimate the value of the uncertainty from 
the discrepancy between the pattern speed that we esti- 
mate and that constrained by (Dehnen 1999). 

To find out if stars in orbits associated with the X\{2) 
family can reach the solar neighborhood, we also need to 
estimate the size of our action variable for the u-anomaly. 
If we assume that u ~ then the epicyclic amplitude is 
set by the condition that we are at r . In other words 
the maximum value in r — r c is approximately rgdr. In 
terms of the coordinate system of our Poincare maps, a 
star which reaches the solar neighborhood has s ~ dr. 

In Figure 2, for dr ~ —0.10 the quasiperiodic regions 
are of size s ~ 0.15 which exceeds the value of dr. This 
is large enough that stars in these orbits would reach the 
solar neighborhood. For dr < 0.093 no Xi(2) orbits ex- 
ist and no stellar orbits aligned perpendicular to the bar 
will reach the solar neighborhood. Because each value of 
dr is associated with a v value in the solar neighborhood 
velocity distribution this implies that above a particular 
value of v no orbits can be in the quasiperiodic orbit region 
surrounding the x\(2) periodic orbits. Our model predicts 
that the w-anomaly should have a sharp boundary in v in 
the solar neighborhood velocity distribution, as observed. 
A model lacking spiral structure would predict this sharp 
edge as well since it is determined by the value of S at 
which the resonance bifurcates (see Figure 1). 

In Figure 2 we see that past a certain value of dr, corre- 
sponding to a value of v, the quasiperiodic orbits are more 
likely to be oscillating with the spiral structure than the 
bar. This implies that there should be an upper boundary 
in v to the Hercules stream or u-anomaly. 

We see in Figure 2 that the quasiperiodic orbits ori- 
ented perpendicular to the bar are bounded by chaotic re- 
gions. This implies that there is a limit in the extent of the 
epicyclic amplitude of stars associated with this quasiperi- 
odic region. In the solar neighborhood, this implies that 
there would be a maximum value of u for orbits in the 
■u-anomaly. Since stars in the chaotic regions can rapidly 
achieve very different epicyclic amplitudes, the boundaries 
of the w-anomaly are probably set by the extent of the 
quasiperiodic orbits associated with the x\(2) orbits. 

3.3. Changing the spiral pattern speed 

According to Contopoulos (1988); Patsis & Kaufmann 
(1999), 2-armed spiral structure should only extend be- 
tween its 2:1 and 4:1 ILRs. These works showed that 
2-armed spiral structure was not supported by the stel- 
lar orbits past the 4:1 ILR. At first we consider pattern 
speeds that are faster, and so with the solar neighborhood 
nearer the 2:1 ILR. The result of integrating Eqn.(24) with 



a faster pattern, f2 s = 0.4, in units of the solar angular ro- 
tation rate f2o is shown in Figure 3. This pattern speed 
would be consistent with a 2-armed spiral pattern near its 
2:1 ILR. 

In Figure 3 we see large scale chaotic regions that com- 
pletely disrupt the spiral structure and the orbits associ- 
ated with the bar OLR. Because v is smaller in this inte- 
gration than that shown in Figure 2, the resonances are 
closer together and so more fully overlapped. The scale 
of the chaotic regions can be reduced by decreasing the 
strength of the bar and spiral perturbations or the values 
of e, [3. However, decreasing both e and (3 by a factor of 
two does not restore the quasiperiodic island about the 
x\(2) family of periodic orbits. The scale of the chaotic 
zone is more sensitive to the extent the resonances overlap 
or to v which is set by the spiral pattern speed, than to 
the size of the perturbations. Since we have adopted real- 
istic values for the bar and spiral arm strength we can use 
this sensitivity to place an approximate limit on the pat- 
tern speeds of spiral structure that passes through the bar 
OLR. We find that the £i(2) are completely disrupted for 
Cl s < 0.45 implying that 2-armed spiral structure passing 
through the OLR is likely to be rotating faster than this. 

In Figure 4 we have integrated a spiral pattern at a 
faster pattern speed of fl s = 0.85 for m = 4. This would 
correspond to an 4-armed spiral pattern. The situation is 
similar to that seen in Figure 2. The OLR bar disrupts 
the spiral pattern on the outer side of the resonance. The 
x\(2) orbit family appears abruptly at dr ~ 0.1, and the 
quasiperiodic island oscillates with the bar until dr ~ 0.12 
where it more strongly supports the spiral structure. 

If we reduce the pattern speed of 4-armed spiral struc- 
ture to fi s < 0.75 phase space looks like that shown in 
Figure 3. Both the x\(2) orbit family and quasiperiodic 
orbits that support the spiral structure are destroyed. This 
implies that 4-armed spiral structure faster than f2 s = 0.75 
cannot pass through the bar's OLR. 

4. SUMMARY AND DISCUSSION 

In this paper we consider the dynamics of stars that 
are affected by perturbations from both spiral structure 
and the Milky Way bar. We construct a simple one- 
dimensional Hamiltonian model for the strongest reso- 
nances in the epicyclic action and angle variables. 

This Hamiltonian is time dependent because the bar and 
spiral structure in general have different pattern speeds, 
and resembles a forced pendulum. To address the issue of 
orbit stability and characterize the nature of phase space 
in this system we numerically integrate this Hamiltonian. 
By plotting points only at the period which separates the 
two resonant perturbations we construct Poincare maps 
which illustrate where in phase space there are area filling 
or chaotic orbits instead of quasiperiodic orbits. 

The Liapunov time in the chaotic regions is of order the 
period which separates two resonances. Over much of the 
resonance this timescale is of order a few times the rotation 
period. Consequently stars which are located in chaotic re- 
gions can move across the region in a few rotation periods. 
This implies that they can achieve large epicyclic ampli- 
tudes and will not remain in any coherent small region in 
phase space. Streams seen in the local velocity distribu- 
tion are unlikely to be located in large chaotic regions. In 
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contrast, particles which are located in quasiperiodic re- 
gions should maintain their epicyclic amplitudes for much 
longer times. Our model predicts that the chaotic regions 
can form large bands. This is a qualitatively different pic- 
ture than the narrow regions spanning a range of Liapunov 
times illustrated by Fux (2001) who considered solely per- 
turbations from the bar. 

We find that larger spiral pattern speeds Sl s cause less 
overlap in the two resonances and so smaller chaotic zones. 
Spiral structure at a pattern speed which puts the solar 
circle is near the 2:1 ILR produces extremely large bands 
of chaos near the 2:1 OLR with the bar. This disrupts 
the spiral pattern and destroys orbits perpendicular to the 
bar near the OLR. It is unlikely that 2-armed spiral struc- 
ture at patterns faster than Q s = 0.45 extends into the 
OLR. Likewise, it is unlikely that 4-armed spiral faster 
than Q s = 0.85 extends into the OLR. It is possible that 
the OLR forms a boundary between spiral structures at 
different pattern speeds. 

2-armed and 4-armed spiral structure at slower pattern 
speeds are much less disruptive. For a 2-armed pattern 
with £l s = 0.6, chaotic regions appear near the OLR which 
disrupt the spiral pattern only across a narrow region. Or- 
bits oriented perpendicular to the bar appear at asymmet- 
ric drift velocities and epicyclic amplitudes consistent with 
existence of quasiperiodic or stable orbits in the local ve- 
locity distribution at the Hercules stream or w-anomaly. 
The boundaries of the it-anomaly are set by the extent 
of the quasiperiodic orbits associated with the family of 
orbits perpendicular to the bar. 

In this paper we have considered the combined prob- 
lem of a resonance from spiral structure and from a 
Galactic bar. Notably the solar neighborhood contains 



a wealth of structure at smaller velocities than the u- 
anomaly (Dehnen 1998). Much of this structure contains 
old stars and so is likely to be caused by spiral structure. 
Since there is more than one clump in the distribution 
it is likely that the solar neighborhood is influenced by 
spiral structure at more than one pattern speed. Hence 
we expect complicated phenomena such as that discussed 
here. In future work we will concentrate on the interplay 
between different modes of spiral structure. 

4.1. Caveats 

Our simple Hamiltonian model neglects the coupling be- 
tween ii , the epicyclic action variable and I2 which is con- 
jugate to the azimuthal angle and related to the angular 
momentum. By neglecting this coupling, we have assumed 
that the particle angular momentum is nearly conserved. 
Notably in Eqn. 1 the coefficient b is of the same order 
as a, implying that neglecting the b term is only a good 
approximation when I2 is small. 

Nevertheless taking the limit of small I2 is reasonable 
first approach and we feel that our simplified model is 
likely to qualitatively illustrate what happens when multi- 
ple resonances are present. We expect that by considering 
the entire 2-dimensional system other phenomena may be 
exhibited, for example particles integrated numerically in 
a system with both spiral structure and a bar can drift in 
both angular momentum and energy. 

This work could not have been carried out without 
helpful discussions with Larry Heifer and Don Garnett. 
A. C. Q. gratefully thanks the Technion for hospitality 
and support during the fall of 2001. 
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Fig. 2. — Poincarc maps made by integrating the Hamiltonian of 
Eqn(24) with m = 2 and with a time step of ^ZL. The panels are 
at different values of dr. The top left panel has the smallest value 
of dr corresponding to a location near the solar neighborhood. The 
lower and righter most panel has the most negative value of dr cor- 
responding to a radius closer to the Galactic center. All panels have 
the same values of e, describing the strength of spiral structure, fi s 
the pattern speed of the spiral structure in units of the solar neigh- 
borhood angular rotation rate (Co), /3, describing the strength of 
the bar, and fif,i the pattern speed of the bar in units of Qq. Phase 
space, as illustrated by the structure in these maps, has two types of 
orbits; curved linear structures, corresponding quasiperiodic orbits, 
and area filling orbits corresponding to chaotic regions. The radial 
distance on these plots is the same as the radial epicyclic amplitude 
in units of ro, or the distance the orbit reaches from r c . A dot is 
placed at the origin for reference. The spiral pattern speed consid- 
ered here would result from a 2-armed spiral pattern with the solar 
neighborhood just within the 4:1 ILR. 



Fig. 3. — Similar to Figure 2 but we have changed the pattern speed 
of the spiral structure to Q s = 0.4. This pattern speed could result 
from a 2-armcd spiral pattern with the solar neighborhood fairly 
near the 2:1 ILR. Since v is smaller the resonances are more fully 
overlapped than is the case for Figure 2. Both the spiral struc- 
ture and resonant orbits at the OLR are disrupted. A resonant 
island of £i(2) orbits do not exist for Q s < 0.45. The stability of 
this orbit family is more sensitive on v or the spiral pattern pat- 
tern speed than the strength of the spiral and bar perturbations. 
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Fig. 4. — Similar to Figure 2 but we have changed the pattern speed of the spiral structure to Q s = 0.85, with the solar neighborhood 
outside the 4:1 ILR. For this simulation m = 4. This system would correspond to a four armed spiral pattern. Phase space as illustrated by 
this map is similar to that shown in Figure 2. 



